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ABSTRACT 

Using samples of structures identified by a multi-scale decomposition from 
0^ I numerical simulation, we analyze the scale-dependence of the virialization of clusters. 

j^ I We find that beyond the scale of full virialization there exists a radius range over which 

clusters are quasi- virialized, i.e. while the internal structure of an individual cluster is 
fj ■ at substantial departure from dynamical relaxation, some statistical properties of the 

^^ ■ multi-scale identified clusters are approximately the same as those for the virialized 

QQ . systems. The dynamical reason of the existence of quasi- virialization is that some of 

CN , the scaling properties of dynamically relaxed systems of cosmic gravitational clustering 

approximately hold beyond the full virialization regime. This scaling can also be seen 
^ I from a semi-analytic calculation of the mass functions of collapsed and uncoUapsed 

pg ' halos in the Press-Schechter formalism. 

^O ' The "individual-statistical" duality of the quasi- virialization provides an explanation 

of the observed puzzle that the total masses of clusters derived from virial theorem 



o 



0^ . are statistically the same as the gravitational lensing determined masses, in spite of 

the presence of irregular configuration and substructures in individual clusters. It also 






% 



explains the tight correlation between the velocity dispersion of optical galaxies and 
the temperature of X-ray emitting gas. Consequently, the virial mass estimators based 



O 

5-H ' on the assumptions of isothermal and hydrostatic model are statistically applicable 

^ ■ to scales on which the clusters are quasi- virialized. In the quasi-virialization regime, 

the temperature functions of clusters also show scaling. This feature is a useful 

discriminator among models. As a preliminary comparison with observation, the 

discriminator yields a favor of the models of LCDM and OCDM. 

Subject headings: cosmology: theory — galaxies: clustering — large scale structure of 
universe — methods: numerical 



1. Introduction 

Clusters of galaxies are playing an important role in the study of large scale structures of the 
universe. In particular, the observed dynamical properties of clusters set stringent constrains on 
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models of structure formation as well as on cosmological parameters. The popular inflation plus 
the cold dark matter (CDM) cosmologies are found to be able to match the existing observations 
of clusters if the cosmological parameters such as the mass density fig, the cosmological constant 
JIa and the mass fluctuation amplitude are properly selected (e.g. Bahcall & Cen 1992; Jing et 
al. 1993, 1995; Jing & Fang 1994; Eke et al. 1996; Viana & Liddle 1996; Carlberg et al. 1997; 
BahcaU, Fan & Cen 1997; Xu, Fang & Wu 1998; Xu, Fang & Deng 1999). 

Yet, some fundamental properties of clusters are not directly measurable, among which the 
conventional determination of gravitating masses of clusters is related to the dynamical state of 
the clusters, i.e. the virialization of both baryonic and nonbaryonic matter. For instance, the 
so-called dynamical mass estimator relies on the hypothesis that the optical galaxies and/or the 
hot intracluster gas are the tracers of the underlying gravitational potential. Nevertheless, the 
virialization assumption is often challenged by the existence of substructures (e.g. Richstone, 
Loeb & Turner 1992; Jing et al. 1995; Crone, Evrard & Richstone 1996; Patrick 1999; Mathiesen, 
Evrard & Mohr 1999). Substructures in both optical and X-ray maps of clusters suggest that 
clusters may still be in the process of formation (e.g. Stein 1997; Kriessler & Beers 1997; Solanes 
et al. 1999). Spatially-resolved measurements of gas temperature in some clusters illustrate 
the complex two-dimensional patterns including the asymmetric variation and the significant 
decline with radius (Henry & Briel 1995; Henriksen &: White 1996; Henriksen & Markevitch 
1996; Markevitch 1996), indicative of strong substructure merging. The unidentified or neglected 
substructures in clusters may lead to an overestimate of galaxy velocity dispersion, and therefore 
the virial mass (Small et al. 1997). This implies that the hydrostatic equilibrium hypothesis may 
be inappropriate to the dynamical state of optical galaxies and intracluster gas for some clusters. 

On the other hand, a statistical comparison of the cluster masses determined from optical 
galaxies. X-ray observations and gravitational lensing shows that there is an excellent agreement 
among different mass estimates on scales much greater than the core radii. The discrepancy of 
these mass estimators appears only inside the central regions of clusters where the local dynamical 
activities of galaxies and cooling flows may become dominant (Wu & Fang 1997; Allen 1998; 
Wu et al. 1998). Since gravitational lensing reveals the cluster gravitating mass regardless 
of the matter compositions and their dynamical status, the consistency between these three 
mass estimators indicates that clusters may be regarded as dynamically-relaxed systems as a 
whole. Further support to this argument comes from the study of the correlations between the 
optically determined velocity dispersion of cluster galaxies, the X-ray measured luminosity and 
the temperature. Namely, although many X-ray clusters show irregular configurations and are 
substructured, there are tight correlations between the velocity dispersion of optical galaxies, the 
temperature and the luminosity of the intracluster gas (Evrard, Metzler, & Navarro 1996; Mohr, 
k Evrard 1997; Wu, Fang k Xu 1998; Wu, Xue & Fang 1999; Mohr, Mathiesen, & Evrard 1999). 

In a word, the observational evidences for whether galaxy clusters are virialized seem to be 
contradictory. In this paper, we intend to highlight this contradiction. We shall show that beyond 
the scale of full virialization there exists a radius range in which clusters are quasi- virialized. 



i.e. the internal structure of individual cluster exhibits a substantial departure from dynamical 
relaxation, while many statistical properties of clusters are approximately the same as those 
of the virialized system. For instance, the virial relation a^ oc M^'^, where ay is the velocity 
dispersion and M is the mass of a system within a given radius, is a good approximation on 
scales larger than virial radius, even though many clusters are irregular in shape and show 
substructures on these scales. Actually, such a "duality" is an indicator of the transition from 
pre-virialization to virialization. This transition can roughly be described by three stages: a) full 
virialization, in which both the individual objects and their statistical properties are virialized; b) 
quasi-virialization, in which the ensemble statistical features are similar to those of the virialized 
systems, but the scattering from the averages is still significant; c) pre-virialization, in which 
both individual objects and their ensemble statistical features are significantly different from 
virialization. 

In hierarchical clustering scenario, the larger the scale, the less the virialization. Thus, the 
dynamical evolution from phase a) to b) or from phase b) to c) can be revealed by studying the 
dynamical properties of the objects on the scale of clusters and beyond. In particular, one can 
address the question: On what scale will the transition of phase a) to b) or phase b) to c) occur? 
We are interested in this issue also because the virialization evolution is model-dependent. It is 
hoped that the scale-dependence of virialization is useful for discrimination of various cosmological 
models. 

The present paper is organized as follows. In §2, we will describe our samples of cluster 
halos which are multi-scale identified from N-body simulations under three different cosmological 
models: 1) the standard cold dark matter (SCDM), 2) low density, flat CDM model with a 
non-zero A (LCDM), and 3) the open CDM (OCDM) model. The Press-Schechter (PS) formalism 
is also introduced so that the simulation results can be confirmed by a semi-analytical approach. 
We then investigate in §3 various statistical indicators of virialization, including density contrast, 
configuration, substructure and velocity dispersion. The emphasis is to find the scale-dependence 
of these features, and determine the scale range of virialization and quasi-virialization. In §4 
we will discuss the applicable range of the virial mass estimator as a statistical measure for 
quasi-virialized systems. In §5 the temperature function of clusters will be constructed and used 
for model discrimination. Finally, we will briefly summarize our conclusions in §6. 



2. Simulated samples of r^rclusters 

2.1. Samples 

The samples to be used for analyzing the scale-dependence of virialization are produced by 
N-body simulations with the P^M code developed by Y.P.Jing (Jing & Fang 1994; Jing et al. 
1995). We employed the following cosmological models: 1) SCDM, 2) LCDM, and 3) OCDM. The 
cosmological parameters (J^m, ^A, ^5 o's) are taken to be (1.0,0.0,0.5,0.62), (0.3,0.7,0.75,1.0), and 



(0.3,0.0,0.75,1.0) for the SCDM, LCDM and OCDM, respectively. These parameters provide a 
consistent description for many observational properties of the universe, especially the abundance 
of clusters (e.g. Jing k Fang 1994; Bahcall, Fan & Cen 1997). 

Other parameters in our simulations, i.e. box size L, number of particles Np and the effective 
force resolution r], are chosen to be {L,Np,7]) = {3Wh~^ Mpc, 64^,0.24/i~^ Mpc). We have run 8 
realizations for each model. A particle has mass of 3.14 x lO^^f^M h~^ Mq, which is small enough 
to resolve reliably the rich clusters of M > 3.0 x 10^^/i~^ Mq. We have also performed simulations 
with a COBE normalization ug = 1 for SCDM. Three realizations were made. 



2.2. Identification of rcrdusters 

There are two ways of studying the scale-dependence of dynamical status of clusters. 
One is to identify clusters with a given size, like Abell radius rab = 1.5 h^^ Mpc, and then 
study the dependence of dynamical properties of these clusters on scales less or greater than 
Tab- Another is to identify clusters by multi-scale resolution analysis, and produce samples of 
rc/-clusters, i.e. clusters identified with scales rd, and then, study the rc^-dependencies, and also 
the radius-dependence of rc^-clusters. 

If all objects on all scales rd are virialized, rcz-clusters should show "clouds in clouds" 
morphology, i.e. all smaller rd clusters locate in the centers of larger rd clusters. In this case, the 
multi-scale identification will actually produce the same samples as the Abell radius identification. 
The Tafe-identification and the multi-scale identification would be essentially equivalent. However, 
in the evolutionary stage of pre-virialization, the two identifications are different from each other. 
The multi-scale identification is more effective to study problems of scale-dependence of clusters 
and beyond (Xu, Fang & Wu 1998; Xu, Fang & Deng 1999; Fang & Xu 1999). 

The detailed procedure of the multi-scale identifications based on the discrete wavelet 
transform (DWT) has been presented in the papers of Xu, Fang & Wu (1998) and Xu, Fang & 
Deng (1999). The technique of applying the DWT for large scale structure study is reviewed in 
Fang & Thews (1998). Briefly, we first describe the distribution of the particles by a 3-D matrix, 
and then do the fast 3-D Daubechies 4 DWT and the reversed transformations to find the wavelet 
function coefficients (WFCs) and scaling function coefficients (SFCs) on various scales. For each 
scale, the cells with SFCs larger than those of the random sample by a given statistical significance, 
say 3a, are picked up as the halos of cluster candidates. Around each of the candidates, we place 
a 6^ grid with the size of cluster diameter and search for the accurate center. The cluster center is 
taken as the position where the largest mass is surrounded. The mass M of a cluster is measured 
by counting the particles within a sphere of radius rd- It is enough to take rd to be numbers of 
bi-fold, i.e. 0.75, 1.5, 3, 6, 12 and 24 /i~^Mpc. Since the rd sphere does not depend on the shape 
of the wavelet functions, this identification is independent of the selection of wavelet function. 

With the identified rc/-clusters, one can find the mass function (MF) or abundance n(> M, rd), 



which is the number density of rc/-clusters with mass larger than M. As an example the MFs of 
rc;-clusters of the SCDM model are ploted in Fig. 1. It shows scaling behavior. For other models 
the MFs show similar behaviors (Fang & Xu 1999). The abundance n(> M,rab) of clusters on 
the Abell radius given by the DWT identification is found to be the same as those given by usual 
friends-of-friends (FOF) identification (Xu, Fang & Wu 1998; Xu, Fang & Deng 1999). 

Since each halo or cluster is characterized by two parameters M and rd, it is inconvenient 
to define richness by its mass alone. Instead, the richness can be defined by the abundance, 
n(> M,rci), or by the mean separation of neighbor clusters 

'^=n(>M,r,0^/3- (1) 

The richness of ri-clusters with masses > Mi is equal to the richness of r2- clusters with masses 
> M2 if their d's are the same. For instance, for the top 25 clusters on various scale rd in the 
simulation box 310'^ h^""' Mpc^, their richness is d = 310/25^'^ ~ 106 h^^ Mpc. For a given r^/, 
this richness definition has its usual meaning, i.e. the larger the M, the less the n(> M,rd) and 
the higher the richness. 



2.3. The Press-Schechter formalism for rcrdusters 

The MFs given by the N-body simulation can be explained by the hierarchical clustering of 
cosmic structures. In the models of hierarchical clustering with Gaussian initial perturbations, the 
comoving number density of collapsed halos at z = in the mass range M to M + dM is given by 
(Press & Schechter 1974) 

where 6c ~ 1.69 almost independent of cosmologies. a{M) is the linear theory rms mass density 
fluctuation in spheres of mass M at redshift z = within a top-hat window of radius R, and is 
determined by the initial density spectrum P{k) and normalization factor ag = A{8h^^Mpc,0). 

From eq.(2), the cumulative number density of the collapsed halos with mass greater than M 
is given by 

n(> M) = / N{M)dM. (3) 

Jm 

Equations (2), (3) have been found to be good descriptions of the mass functions of collapsed 

clusters identified by FOF method from N-body simulation (e.g. Mo, Jing & White 1997). Figure 

1 shows that PS predications fit very well with the mass function (MF) of clusters on the Abell 

radius given by the DWT identification. 

As pointed out by Mo & White (1995) and Bi & Fang (1996; 1997), the PS formalism can 
be employed to approximately describe not only collapsed halos, but also uncoUapsed regions. 



i.e. both collapsed halos and uncollapsed regions can be viewed as the sum of various individual 
top-hat spheres. An uncollapsed region (or "quasi- virialized halo" in terminology of this paper) 
corresponds to a region in which the initial linear fluctuation is less than the threshold 1.686 at 
the redshift considered. As a consequence, one can consider the initial density field as a system 
consisting of many spheres, each of which has initial radius R and density contrasts 5q, and mass 
is equal to M = Ivri^^po- 

Let's consider an arbitrary spherical volume of radius r^i at redshift z = 0. Matter in this 
volume can come from various initial spheres with R and 60. For a given r^i, one can find the 
relationship between R and Sq by considering the nonlinear evolution of a massive shell under 
spherical symmetry (Padmanabhan 1993; Mo & White 1995; Bi & Fang 1996). Then, we have the 
relation 

|=/(^). (4) 

It means that a spherical region of 6th ^-nd R will evolve into r^i at z = 0, which includes both 
collapsed and uncollapsed objects. The function f{x) is given in Bi & Fang (1996). 

If we identify clusters with radius Vd in an arbitrary spatial domain Ax, the fraction of the 
total mass poAx contained in clusters of mass > M corresponds to the fraction of mass which 
have initial fluctuations larger than a given 6th- For a Gaussian field, this fraction is given by 

/■°° 1 f 6^ \ 

F(> M, r,;) = / -^ exp --^ d6, (5) 

where a\j is the variance of density perturbation on scale R. Because a spherical region of 6th and 
R will evolve into r^i at z = 0, initial spherical regions oi 6 > 6th and R will evolve into radii less 
than Tc/ at z = 0. Therefore, the comoving number density of the collapsed and uncollapsed PS 
spheres with mass in M ^ M + dM, which are spheres identified with radius rd at z = 0, is given 
by 
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dM (6) 



2 po -^/ dlna dln6 
^M2^^ ' ^dlnM ~ dlnM'' 

where v = 6{R,rci)/(T(R) denotes the relative height of the initial perturbation (evaluated at 
redshift z = 0). The number 2 in the equation is the notorious ad hoc multiplication factor of PS 
formalism (see discussion in Bond et al. 1991). 

The cumulative number density, n(> M,rci), of collapsed and uncollapsed rc^-halos with mass 
greater than M should be 

n{>M,rd) = / N{M,rd)dM (7) 
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These MFs are also plotted in Fig. 1. It shows that the PS formalism basically explains the 
mass functions given by N-body simulation, especially the scaling behaviors. Meanwhile, Fig. 1 
shows that the PS calculations systematically underpredict the MFs on scale r^ > 3 h~^ Mpc. 
This is most likely caused by the hypothesis of spherical symmetry of the dynamic evolution of 
uncollapsed regions. A similar discrepancy has also been found in the two-point correlation of 
halos in underdense (or uncollapsed) regions, i.e. the PS formalism significantly underpredicts the 
halo-halo correlation when small halos are in low-density environments (Jing 1998). Nevertheless, 
the PS result shows the scale-dependence of statistical features given by the PS formalism is 
the same as that of N-body simulation. This agreement is very useful to reveal the status of 
quasi-virialized (or uncollapsed) regions. 



3. Virialization and quasi-virialization 

3.1. Density contrasts of rc/-clusters 

A popular indicator of virialization of clusters is the density contrast. For spherical collapse 
in a standard flat universe, the virial radius is usually measured as the radius within which the 
mean density contrast is ~ 178, i.e. Pvir/Pb = 178. For open universe, the density contrast within 
virial radius is in the range of ~ IOO/JIm to 200/^1^ (Lacey & Cole 1993). For models of SCDM, 
LCDM and OCDM, the density contrasts Pvir/pb a-re equal to 178, 335, 402, respectively. 

Observationally, the spatial range of virialization is characterized by r2oo, the radius where 
the average interior number density of galaxies is 200 times higher than the average. Thus, 
the deviation from virialization of a rc^-cluster can be quantified by the mean mass contrast, 
(AM)/(M), where (M) is the mean mass within radius rd, and (AM) the mean mass excess of 
the rc^-clusters. Fig. 2 plots the mass contrasts of rc/-clusters with richness d = 30, 50 and 100 
h-^ Mpc for the SCDM, LCDM and OCDM models, respectively. The lines in LCDM almost 
overlap with those for OCDM. As expected, the more massive (i.e. the larger d) the clusters, the 
closer to virialization. For the SCDM clusters with richness d = 30, 50, and 100 h~^ Mpc, the 
virialized radii r^r at z = are 1.09, 1.33, and 1.60 /i~^ Mpc, respectively. They are 1.03, 1.33, 
1.82 h~^ Mpc for LCDM, and 0.96, 1.27, 1.70 /i^^ Mpc for OCDM. 

Therefore, except for the richest clusters, the full virialization generally is realized only on 
scales less than rat = 1-5 h~^ Mpc. This is consistent with the observed result that r2oo is usually 
smaller than the Abell radius rat- Namely, except for the core regions of clusters, Abell clusters as 
a whole are not yet fully virialized. 

In the case of positive bias, r2oo given by galaxies is larger than that of dark matter, and 
therefore, the real virialized range will be even smaller than that given by galaxy observations. As 
a result, in the scale range of equal to and larger than the Abell radius, halos for rc^-clusters are 
mostly in the stages of quasi- or even pre-virialization. 



3.2. Configurations of rcrdusters 

The spatial configuration of clusters is also a useful indicator of virialization. As discussed 
above, if the halos of clusters are fully virialized, their spatial configurations have to be regular 
"clouds in clouds", i.e. the smaller rc/-clusters are located in the centers of larger rc^-clusters. 

Fig. 3 plots a projected distribution of rc^-clusters identified from one realization of model 
OCDM within box 310 h^^ Mpc^, which illustrates the top 25 massive clusters for each r^i (the 
richness \s d = 106 h"^ Mpc). The scales r^i are taken to be 0.75, 1.5, 3, 6, 12 and 24 /i^^Mpc, 
respectively. Fig. 3 shows a variety of configurations of the halos. Apparently, some halos show 
"clouds in clouds" morphology, while some do not. 

According to the scale of r2oo in Fig. 2, it is expected that, for such a high richness, clusters 
of Tci < 2 h~^ Mpc have to be symmetrical and concentric "clouds in clouds" structures. Fig. 3, 
indeed, shows that 21 of the 25 r^; = 0.75 and 1.5 /i~^ Mpc are concentric. These r^i = 0.75-1.5 
"clouds in clouds" objects are fully virialized. 

On the other hand, there are seven Vd = 3 h^^ Mpc clusters which do not contain either 
Vci = 0.75 h~^ Mpc or 1.5 h~^ Mpc clusters in their centers. These are massive clusters on scale 
3 h^^ Mpc, but without fully virialized cores. On scales r^ > 3 h^^ Mpc, more rc^-clusters are 
asymmetric, irregular, and without small-scale high density peaks within them. These features 
show that identifying structures with a given scale like Vat may miss massive objects which do not 
contain high density peaks on scale rat- Only in the range r^i < T20O) the multi-scale identifications 
find the same sample as the single-scale identification. 

In terms of mass and length scale, rd = 24 h^^ Mpc clusters actually are superclusters. 
Observationally, supercluster is defined as cluster of rich clusters (top Tab-clusters). In Fig. 3, one 
can see that some rah-clusters are located together to form larger scale clumps or filaments. They 
are obviously superclusters. Fig. 3 shows, however, that some of the top 25 r^i = 24 h~^ Mpc halos 
contain nothing of r^i < 24 h^^ Mpc massive halos. This is, although the mass and length scale 
of these objects are the same as superclusters, they do not contain any of the top 25 clusters on 
smaller scales. This result indicates that the identification by "cluster of rich clusters" probably 
is incomplete. Recently, it has been found that the dynamical properties of superclusters are not 
only determined by rich clusters, but also by numerous galaxies in the intra-cluster regions (Small 
et al. 1998). Therefore, to study the dynamics on scales beyond clusters, one should directly 
identify structures from matter or galaxy distributions, instead of using only rich clusters. 



3.3. Substructures and potential minimums of rc/-clusters 

Similar to the configuration, the number of substructures in a cluster is also an effective 
measure of the deviation from virialization. The asymmetry of the configuration of clusters 
actually is caused by their substructures. 



Considering that the time scale required for a sound wave in the intracluster gas to cross 
a cluster is shorter than the dynamical time of the cluster, it is reasonable to assume that the 
intracluster gas is in hydrostatic equilibrium with local gravitational potential, and therefore. 
X-ray images map the cluster potential locally (Jing et al. 1995). Thus, one can identify 
substructure as potential minimum of clusters. 

We search for the potential minimums by the same algorithm of Jing et al. (1995). Briefly, we 
first place a mesh with 0.24 /i^^Mpc resolution around each center of identified rc^-cluster halos. 
The size of this mesh is taken to be 7.5, 7.5, 7.5, 15 and 30 h^^Mpc for clusters of r^i = 0.75, 
1.5, 3, 6 and 12 /i^^Mpc, corresponding to grid points 32^, 32^, 32'^, 64^ and 128'^, respectively. 
We accumulate mass density values for the meshes from particles. The mass of each particle is 
smoothed by a Gaussian kernel 

where s is the smoothing length. The smoothing length Sj of particle i is equal to the local mean 
separation di, calculated by counting five nearest particles around i. Then the density value on 
arbitrary cell j is given by 

Pj = ^miW{rji,Si), (9) 

i 

where Vij is the separation between cell j and particle i of mass rrii. The summation in Eq.(9) is 
made over all particles within radius 3.75 h^"*^ Mpc. With this mass distribution, gravitational 
potential on the grid can be obtained. We then identify a cell to be a potential minimum if its 
potential value is smaller than those of its all 26 neighbors. 

Obviously, if a rc/-cluster contains only one potential minimum, it means no substructure. 
Fig. 4 shows the dependence of average number of substructures per rc^-cluster on the count of 
the particles within the cluster, which is proportional to the cluster mass. As expected, the larger 
rd clusters contain more substructures on average because larger scale objects are less virialized. 

A "surprising" result shown in Fig. 4 is that for r^i > 3 h^^ Mpc clusters, the more massive 
the clusters (i.e. the more the count of particles), the more the substructures. This seems to 
be contradicting with the general conclusion of last section: the more massive the cluster, the 
closer to virialization. Actually, this is due to the threshold used for identifying substructures. In 
massive objects substructures with deep potential valleys are easy to grow up, but not so for less 
massive objects. For clusters with rd = 1.5 /i~^Mpc, containing no substructures on average (as 
shown in Fig. 4) doesn't mean that all r^ = 1.5 /i^^Mpc clusters are fully virialized. Contrarily, 
low richness r^ = 1.5 /i^^Mpc clusters are pre-virialized (Fig. 2). They simply do not have enough 
mass to form substructures above the defined threshold. 

Similar phenomena can be seen from Fig. 5 which shows the probability distributions of the 
number of substructures per rd = 24 h^^Mpc clusters with richness d = 80 h^^ Mpc, i.e. the top 
59 richest clusters in the simulation box 310^ /i^^ Mpc^. From Fig. 5 it is seen that the SCDM 
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clusters generally are less substructured than those of the LCDM and OCDM. It doesn't mean 
that the SCDM clusters are more close to virialization. Instead, the deep potential valleys for 
substructures in the SCDM r^i = 24 h~^Mpc clusters formed later than the LCDM and OCDM. 

Fig. 6 gives the rc;-dependence of the mean number of substructures per rc^-clusters having 
richness d = 80 h~^ Mpc. For all the three cosmological models, the number of substructures per 
rc/-clusters increases with Vd- As Fig. 5, the mean number of substructures in model SCDM is less 
than those of LCDM and OCDM. The substructure number at redshift z = 0.5 is less than that 
of z = 0, too. This is once again to indicate that deep potential valleys have not yet developed at 
z = 0.5. 

Thus, despite substructure is an indicator of the deviation from virialization, there is no simple 
relation between virialization and substructures. A system with rich substructures definitely 
deviates from virialization, but the existence of substructures is not a necessary condition of the 
deviation from virialization. Systems with more deep potential valleys are not always in a lower 
degree of virialization than those with fewer potential valleys. This point is important in studying 
objects on scales larger than Tab- 



3.4. Velocity dispersion 

The most decisive measure of virialization certainly is the velocity dispersion of a cluster. For 
a virialized rc^-cluster we have o"^ oc M(< rd), in which M(< rd) denotes the mass of a cluster 
within radius r^i- The velocity dispersion a^ of rcrduster is calculated from the velocities of all 
CDM particles within radius r^i with respect to the center of the cluster. Fig. 7 shows the three 
dimensional a^ of r^i = 0.75 — 24 h~^ Mpc clusters identified from three realizations of the OCDM 
model at z = 0. 

In the case of full virialization, a^ of the system with a given size rd should be completely 
determined by its masses M(< rd)- Fig. 7 shows that a^ is almost completely determined by rd 
and M(< rd) if clusters are as massive as M(< rd) > 10^^'^, 10^^ and lO^^'^/i^^M©, respectively 
for rd = 0.75, 1.5 and 3 h~^ Mpc. Namely, in these parameter ranges (r^;, M), clusters are fully 
virialized. Out of it, velocity dispersion a^ is no longer completely determined by radius r^/ and 
mass M{< rd)- As shown in Fig. 7, the scattering of velocity dispersion a^ around their mean a^ 
can be as large as 15% (see the vertical error bars of Fig. 7). This scattering is not caused by 
Poisson fluctuations. For rd = 0.75 — 3h~^ Mpc, the scattering increases with the increase of rd 
and the decrease of M(< rd)- Therefore, the scattering is due to the deviation from virialization. 

The mean velocity dispersion a^ of the panels of rd < 6 h~^ Mpc in Fig. 7 is found to be 
proportional to M^/^(< r^/). Fig. 8 plots the result of 7 = (ilog(T„/dlogM(< rd) for the data set 
shown in Fig. 7. Fig. 8 clearly demonstrates that for scales rd < 6 h~^ Mpc and richness d = 
20 - 100 h^^ Mpc, all the three models yield 7 ~ 0.5, which is the value for virialized systems. 
Namely, the ensembly averaged relations between mass and velocity dispersion of rd < 6 h^^ 
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Mpc clusters are the same as those for the viriahzed systems. Even for ensemble of rich r^ = 12 
h^^ Mpc clusters with d > 80 h^^ Mpc, virialization holds approximately true. However, this 
approximation is no longer valid for the rd = 24 h~^ and d < 100 h^^ Mpc clusters. 

Strictly speaking, virial theorem is i^ = W/2, where W and K represent, respectively, the 
potential and kinetic energy of the system. Therefore, it includes two numbers 1.) the index 
7 = 0.5 from K oc W; 2. the coefficient 1/2 from K/W = 1/2. In the above used parameter 
space of r^i and d, both the index and the coefficient are found to be not different from 1/2 by 
10%. One can define "quasi- virialization" to be a dynamically evolutionary stage of rc/-clusters, 
for which ensemble averaged properties approximately agree with those required by virialization, 
while individual clusters may significantly depart from virialization. 

Thus, Figs. 7 and 8 show that the quasi-virialization range for the LCDM and OCDM clusters 
at z = is from radius r^j^ to about 6 h^^ Mpc with richness d > 20 h~^ Mpc, and to about 12 
h-^ Mpc with d > 80 /i^^ Mpc. For the SCDM clusters of r^ ~ 12 h'^ Mpc, only the richest 
clusters {d > 100 h~^ Mpc) are quasi- virialized. The value 7 at redshift z = 0.5 is also plotted 
in Fig. 8. For the LCDM and OCDM, the quasi-virialization is already realized for structures of 
Vci < 6/i~^Mpc at z = 0.5. For the SCDM, the quasi-virialization approximately holds for clusters 
of rd < 6 h^^ Mpc, while rd = 12 h^^ Mpc clusters are completely out of the quasi-virialization 
range. 

The redshift evolution of 7 for models SCDM, LCDM and OCDM is shown in Fig. 9, which 
gives the mean 7 at redshifts z = 0, 0.2, 0.5, 0.8 for samples of top 59 clusters, or richness d = 80 
h~^ Mpc, where rd is taken to be 1.5, 12 and 24 h^^ Mpc. For clusters of rd < 12 h^^ Mpc in the 
LCDM and OCDM models, we have 7 ~ 0.5 in the redshift range of z < 0.8. But in the SCDM 
model, the stage of quasi-virialization has not yet been reached at redshift z > for clusters of 
rd = 12 h'^ Mpc. 



4. The virial mass estimator 

With the above analysis on quasi-virialization, one can now study the applicability of the 
virial mass estimator. Especially, the individual-statistical (or individual- ensemble) duality 
of quasi-virialization is very useful for determining the available range of the virial mass as a 
statistical measure. 

The virial estimator of the total mass of a cluster within radius r is given by (e.g. Bahcall & 
Sarazin 1977; Mathews 1978), 

where pgai is the density of galaxies. For isothermal sphere approximation of the galaxy 
distribution, we have /(r) = 3r^/(r^ -|- r^), where re is the core radius. Basically, the applicability 
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of eq.(lO) relys on the following three assumptions: 1.) The second term in the bracket of eq.(lO), 
i.e. dlna'^/dlnr, is negligible. Namely, o"^ should be independent of r; 2.) For a given r, the 
velocity dispersion o"^ is proportional to M(< r), i.e. 7^, = dlnay/dlnM = 0.5; 3.) The profile of 
mass distribution /(r) doesn't significantly depend on a^ and richness d. 

The assumption 2 has been studied in §3.4. The mean 7 is indeed equal to 0.5 for 
quasi-virialized clusters. We will show below that the assumptions 1 and 3 are also statistically 
available for quasi-virialized systems. 



4.1. Profiles of velocity dispersion 

To test the assumption 1, we calculated the mean velocity dispersion a^ within radius r 
among the 50 most massive clusters (or richness d '^ 85 h^^ Mpc) at scales rd = 0.75, 1.5, 3, 6, 12 
and 24 h~^ Mpc. The results averaged from 8 realizations for OCDM and SCDM are plotted in 
Fig. 10. The LCDM gives similar results as the OCDM. Fig. 10 shows that the profiles of velocity 
dispersion of r^i < 6 h~^ Mpc clusters are rather flat, and the variations of the velocity profiles 
over the radius range 0.5 < r < 10 h~^ Mpc are no more than 20%. Beyond r = 10 h~^ Mpc, the 
velocity of particles is dominated by the Hubble expansion. 

Fig. 10 also shows that a^ is virtually rc/-independent for clusters of r^ < 3 h^^ Mpc. When 
Tci > 3 h^^ Mpc, a^ is lower for larger rd- However, the velocity dispersion of rd = 6 h~^ Mpc 
clusters is only ~ 20% lower than that of rd < 3 h^^ Mpc clusters. Even in the case of rd = 12 
h^^ Mpc, the o"^ is not different from r^; < 3 h~^ Mpc clusters by a factor of 2. 

Despite many top rd = 6 and 12 h^^ Mpc clusters do not contain top rd < 6h^^ Mpc clusters 
as "clouds in clouds" , the mean profile of the velocity dispersion (T„ remains to be flat from 1 to ~ 
10 /i"^ Mpc for clusters of rd < 6 h'^ Mpc and d > 40 /i"^ Mpc or rd < 12 h'^ Mpc and d > 80 
h^^ Mpc. Therefore, the assumption 1 is statistically correct for quasi-virialized clusters. 



4.2. Profile of mass distributions 

To test the assumption 3, we calculated the cross-correlation function between clusters and 
mass particles, icm{T~) = N{r)/Nexp{T) — 1, where N{r) is the number of particles within a spherical 
shell with radius ~ r, and Ngxpir) is the mean number expected from random distribution. 

Fig. 11 gives the cluster-particle cross-correlation functions of clusters with rd = 1.5, 6 and 
24 h^^ Mpc and richness d = 90 h~^ Mpc. It shows that the mean cross correlation functions 
Ccmif) are about the same for all rcrclusters. One can describe the correlation function ^cm as 
Ccm oc r"~'^ (or, /(r) oc r") in the radius range from 1 to 20 h^^ Mpc. The index a and their 
variance of each rd are listed in Table 1. For the LCDM and OCDM, all the indexes of rd < 6 
h~^ Mpc clusters are a ~ 0.5. For the SCDM a is slightly lower than 0.5 at r^/ < 6 h~^ Mpc. 
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Fig. 12 plots £,cm{f) for the r^i = 6 /i^-^Mpc clusters with richness from d = 40 to 120 h^^ 
Mpc in the OCDM. icm{f) shows a systematic increase with d, indicating that the mean mass of 
clusters increases with d. However, the shapes of density profiles for all richness can approximately 
be written as ^cm{f) ot r°~^, or /(r) oc r", and a ~ 0.5 in the radius range 1-20 /i~^Mpc. That 
is, the shape of mass profile is independent of the richness. 

Moreover, for a given r^i such as rci= 6 h^^ Mpc, the mean velocity dispersion a increases 
with mass or d (Fig. 7). Therefore, the cluster-particle cross-correlation functions E,cm{f) with 
different d in Fig. 12 actually is for different a. Therefore, Fig. 12 also shows that the shape of 
mass profile also does not depend on the mean velocity dispersion. 

Thus, one can conclude that the mass profile /(r) oc r" with a ~ 0.5 is weakly dependent on 
Tc/, richness and velocity dispersion. One may safely employ the assumption 3 at least on scales 
< 6 /i~^ Mpc, i.e. in the range of quasi- virialization. Actually, the independence of the cluster 
mass profiles on the richness and velocity dispersion has been found by Jing et al. (1995). 

It should be pointed out that the radius range shown in Figs. 11 and 12 are generally larger 
or much larger than r2oo- That is, we study only on scales much larger than the cores of clusters, 
and therefore, the core radius is less important. 



4.3. The goodness of the virial mass estimator 

Accordingly, mass estimator eq.(lO) is good as a statistical measurement for clusters on radius 
larger than rmr or r200) ^-nd within the range of quasi- virialization. The goodness of this statistical 
measurement can be seen from Fig. 7, in which the vertical error bars show the deviation of log a^ 
from their virialization average. The error bars Alogcj^ are of the order of ~ 0.15. Thus, Mj, given 
by estimator (10) should have an uncertainty of about AM^/Mj, ~ 0.30. This result provides an 
explanation for the fact that the cluster virial masses determined from optical galaxies statistically 
are in a good agreement with gravitational lensing-derived masses on scales much greater than 
the core radii. For instance, based on a sample of lensing and optical clusters, the correlations 
between the virial mass M^ and weak lensing-derived mass Mje„ of clusters in the scale range of 
0.25 - 1.5 h-^ Mpc are found to be (Wu k Fang 1997; Wu et al. 1998) 

Mien = (0.63 ± 0.35)M^. (11) 

Since the samples used for eq. (11) are not statistically as uniform as simulation sample, 
one cannot take a detailed comparison of observed results with simulations. Nevertheless, the 
individual-ensemble duality of quasi- virialization gives a reasonable explanation for the eq.(ll) 
which used the sample consisting of pre-dynamically-relaxed clusters. 

Similar to eq.(ll), the mass Mx estimated from hydro dynamical equilibrium model of X-ray 
gas is found to be, on average, in agreement with the gravitational lensing-derived cluster masses 
Miens- The correlation between Mx and Mien in the scale range of 0.25 — 1.5 h^^ Mpc is (Wu & 
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Fang 1997; Wu et al. 1998) 

Mzen = (0.97±0.44)Mx. (12) 

This result also shows the feature of quasi-virialization: the hydro dynamical equilibrium 
model-derived Mx is reliable, while individual clusters may deviate from virialization. 

Moreover, the most recent statistics using the largest sample of 149 clusters gives (Wu, Fang 
& Xu 1998) 

where the parameter fi = 0.59 is the mean molecular weight for intracluster gas and nip is the 
proton mass. Eq.(13) means that the temperature T of X-ray emitting intracluster gas statistically 
can be used as a virial indicator of the underlying gravitational potentials of clusters, and is 
related to the velocity dispersion by virial relation as 

3 1 

-/cT = -fimpa^. (14) 

The current observations have revealed a constant temperature profile within virial radius (Irwin, 
Bergman & Evrard 1999). Since the profile of a^ is flat within radius r < 6 h^^ Mpc (Fig. 10), 
eq.(14) indicates that the averaged temperature of clusters is also approximately constant within 
the range of quasi-virialization. 



5. Temperature function of clusters 

5.1. Temperature functions and dynamical state 

For the virialized and quasi-virialized clusters, we can define the temperature function (TF) 
of clusters, nT{> T, rd), which is the number density of rc/-clusters with temperature larger than 
T. With eq.(14), the TFs can be calculated from the simulated velocity dispersion function (VDF) 
n„(> ay,rci), which is the number density of rcrclusters with velocity dispersion larger than a^. 
Actually, eq.(14) has been widely employed in the calculation of TF from VDF (e.g. Klypin & 
Rhee, G. 1994; Jing & Fang 1994). Yet, as we have emphasized, we are interested in the effect of 
dynamical state on the TFs. 

Fig. 13 plots the differential TFs, n{T,rci) = dn{> T,rci)/dT, for r^i = 0.75 to 12 h'^^Mpc 
clusters in the three dark matter models. In fact, the VDF-TF transfer [eq.(14)] may not be 
correct for clusters of r^ = 12 h^^ Mpc, which is plotted only for comparison. 

If the system is fully virialized on scales r < rd, the TFs should be basically scale- independent 
in this scale range, because virial temperature T is independent on radius, and clusters identified 
by {T,rci) is the same as {T,r'^i). Fig. 13 shows the significant scale-dependence of the TFs, even 
at Vci = 0.75 h^^ Mpc. Therefore, these systems deviate from full virialization. 
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However, Fig. 13 shows that all the TFs of rd < 6 h^^ Mpc have similar shape, i.e. the TFs 
are scaling. This feature is from quasi-virialization. The TFs of r^ > 12/i~^ Mpc clusters do not 
join the scaling of r^ < 6 h~^ Mpc, and therefore, it is not in quasi-virialization state. This is 
expected, as Fig. 7 has shown that rd > 12 h~^ Mpc clusters are different from that of < 6 h~^ 
Mpc. 



5.2. The size and virialization of X-ray clusters 

As an application of simulated TF n{T,rd) (Fig. 13), we analyze the problem of the size 
of X-ray clusters. Observationally, X-ray clusters are described by their X-ray luminosity, flux 
weighted temperature and morphology. Although for some individual clusters we have known their 
X-ray isophotes and therefore the aperture of the X-ray emission, the X-ray luminosity is generally 
published as pseudo-total luminosity which does not contain information about the spatial scale 
of clusters. The Abell radius rab is not used as a condition in X-ray cluster identification. Even 
for Abell clusters, there is also the ambiguity of whether the X-ray luminosity comes from cluster 
halos of larger than radius Tab = 1-5 h~^ Mpc. If clusters are fully virialized, the size uncertainty 
of the TFs is not important since the TFs are size-independent. 

However, in quasi-virialization state, the size-dependence of TFs is substantial. For instance, 
we ploted the observed TF given by Edge et al. (1990) and Henry & Arnaud (1991) in Fig. 13. 
The LCDM and OCDM TFs are in good agreement with the observed data only if rd is in the 
range 0.75 - 1.5 h~^ Mpc. In other words, the LCDM and OCDM can pass the TF test if most 
X-ray clusters used for the observed TF are on the scales of 0.75 and 1.5 h~^ Mpc. On the other 
hand, to fit the SCDM TF with observed data, the size of X-ray clusters should be in the range 
1.5 -3 h^^ Mpc, which is larger than the results of LCDM and OCDM by a factor of 2. 

The role of the size of X-ray clusters can be more clearly illustrated by SCDM model with 
normalization amplitude erg = 1. Fig. 14 shows the TFs of the cg = 1 SCDM model. In this case, 
the model-predicted TF can fit the data only if the size of X-ray clusters is in the range 6 to 12 
h~^ Mpc, which seems very unlikely. 

Thus, one can conclude that for either the model of SCDM or LCDM and OCDM, the 
simulated TFs can fit with observations, if only the size of the X-ray clusters used for constructing 
the TFs is, on average, larger than their virial radius. This is consistent with the fact that the 
configuration of X-ray clusters generally is irregular and substructured. 

Obviously, the scale-dependence of the TFs can be employed for discriminating among models 
if the information of the size of X-ray clusters is available. Currently, we are lacking of the data of 
X-ray cluster sizes. Instead, one can consider the test of TF plus two-point correlation function. 
As it has been shown in Figs. 13 and 14, the observed TFs can set a lower limit to the scale of 
clusters. For instance, the data of Edge et al. (1990) require that rd should be > 0.75 /i~^ for the 
LCDM and OCDM, > 3 h'^ Mpc for the SCDM, and > 6 /j-^ Mpc for the as = 1 SCDM. On 



16 



the other hand, for a given abundance, the larger the scale Vd, the smaller the correlation length 
(Xu, Fang & Deng 1999). At the abundance d > 70 h^^ Mpc, the correlation length is found to be 
ro > 15/i"i Mpc if only r^ < 6 h'^ Mpc in the LCDM and OCDM; and r^ < 3 /i'^ in the SCDM. 

Since the observed correlation length tq of X-ray clusters is rg > 15h~^ Mpc. The LCDM and 
OCDM are able to pass the two constraints to rd simultaneously: 1.) TF requires rd > 0.75 h^^ 
Mpc, and 2.) correlation length requires rd < 6 h~^ Mpc. Yet, the SCDM with ug = 1 is clearly 
in trouble as 1.) TF requires rd > 3 h^^ Mpc, and 2.) correlation length requires at least rd < 3 
h~^ Mpc. 



6. Conclusions 

Using N-body simulated samples in models of the SCDM, LCDM and OCDM, we have 
analyzed the dynamical state of the DWT multi-scale identified clusters. We showed that state 
of quasi-virialization exists in the dynamical evolution of these clusters, i.e. while the internal 
structure of individual clusters significantly differs from dynamical relaxation on scales larger 
than virial radius r^r, some statistical features of the clusters are approximately the same 
as those for the virialized systems. Actually, the virial statistical features, like a-o oc M^'^, 
Ccmi''') Of r^^'^, are scaling relations. The dynamical reason of the existence of quasi-virialization is 
that some of the scaling properties of dynamical relaxed systems of cosmic gravitational clustering 
approximately hold beyond the fully virialization regime. The scaling behavior of mass functions 
of the quasi-virialized clusters can also be repeated by semi-analytic calculation on collapsed and 
uncollapsed halos in Press-Schechter formalism. 

Consequently, in terms of statistical description, the cluster galaxies and X-ray gas can 
be used as virial tracers of cluster potential on scales beyond the full virialization. This result 
provides a good explanation that the virial masses and X-ray masses are basically the same as the 
gravitational lensing determined masses on scales beyond the full virialization. It also explained 
the tight correlation between the velocity dispersion of optical galaxies, and the temperature of 
the clusters. The virial mass estimator based on the assumptions of isothermal and hydrostatic 
model is applicable on scales as large as about rd = 6/i~^Mpc. 

In the quasi-virialization state, the TFs have scaling. It is potentially important for 
model-discrimination. A very preliminary result given by the test of the scale-dependence of the 
TFs showed that the LCDM and OCDM are favored. 

Wen Xu thanks the World Laboratory for a scholarship. This project was done when W.X. 
was visiting at Physics Department, University of Arizona. 
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Table 1. Mean number of a = dlnM/dlnr in the range from 0.5 to 24 h ^Mpc 



rci 0.75 1.5 3. 6. 12. 24. 



OCDM 


0.50±0.06 


0.49±0.06 


0.52±0.05 


0.53±0.05 


0.68±0.04 


0.67±0.06 


LCDM 


0.51±0.06 


0.47±0.06 


0.51±0.05 


0.55±0.04 


0.67±0.04 


0.65±0.06 


SCDM 


0.40±0.08 


0.48±0.06 


0.44±0.06 


0.49±0.07 


0.67±0.06 


0.75±0.07 
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Fig. 1. — The mass function of clusters with rd = 1.5,3, 6, 12, 24/i^-^ Mpc of N-body simulation 
(circle) and of the PS formalism predictions (solid line) for SCDM. For the PS results, the MF of 
r^i = 1.5h~^ Mpc is the number density of collapsed halos, and the MFs of r^ = 3, 6, 12, 24/i~^ 
Mpc are the number density of collapsed and uncollapsed rc;-halos. 

Fig. 2. — Logarithmic mass contrast as a function of r^i for rc/-clusters with richness d = 30, 50 
and 100 h'^Mpc. 

Fig. 3. — A projected distribution of rc^-clusters identified from one realization of OCDM N-body 
simulation within box 310 /i^^ Mpc^. It shows the top 25 of massive clusters for rd = 0.75, 1.5, 3, 
6, 12 and 24 h~^ Mpc. The legends of the symbols are: plus for 0.75, diamond for 1.5, circle for 
3, triangle for 6, square for 12 and giant circle for 24 h~^ Mpc. The linear size of the symbols is 
roughly equal to the corresponding rd, but slightly amplified for smaller r^/ for clarity. 

Fig. 4. — Average number of substructures per rc/-cluster as a function of logarithmic counts of 
particles in each cluster. The data is from one realization of the OCDM model simulation. 

Fig. 5. — Probability distributions of number of substructures per cluster for top 59 richest 
rd = 24/i~^Mpc clusters in a 310 h^^ Mpc simulation box. Results of three cosmological models 
are shown. 

Fig. 6. — The mean number of substructures per cluster as a function of rd for top 59 richest 
clusters in a 310 h~^ Mpc simulation box for models of the SCDM, LCDM and OCDM. Error bars 
are the standard deviation in N. The results of SCDM and OCDM at redshift z = 0.5 are also 
shown, but without error bars for clarity. To avoid overlapping with the line of SCDM at z = 0, 
the line for LCDM at z = 0.5 is not shown. 

Fig. 7. — Three dimensional velocity dispersion a^ vs. mass of clusters of rd = 0.75, 1.5, 3, 6, 
12, and 24 h^^ Mpc. The data are from three realizations of OCDM model at z = 0. The mean 
logarithmic velocity dispersion, logo"^, within each bin of width 0.1 dex in log M(rc/) are shown by 
thick crosses, for which the horizontal bars represent the bin widths, and the vertical error bars are 
the standard deviation of log a^ within each bin. 

Fig. 8. — 7 vs. richness (d) of clusters of rd = 0.75, 1.5, 3, 6, 12, and 24 h^^ Mpc. The results 
for three models SCDM, LCDM and OCDM at redshifts z = and z = 0.5 are given. Error bars 
are only illustrated for rd = 12 h^^ Mpc. Errors for other rd are similar to rd = 12 h^^ Mpc, but 
having not been illustrated for clarity. 

Fig. 9. — The redshift evolution of 7 for the top 59 clusters of rd = 1.5, 12, and 24 h^^ Mpc in 
models SCDM , LCDM and OCDM. Error bars are illustrated only for rd = 12/i"^ Mpc. 

Fig. 10. — The averaged profiles of velocity dispersion within radius r for top 50 massive clusters in 
a 310 h-^Mpc simulation box. Clusters of rd = 0.75, 1.5, 3, 6, 12, and 24 /i^i Mpc of the SCDM 
and OCDM are shown. 
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Fig. 11. — The cross-correlation function between cluster and its particles for clusters with rci= 
1.5, 6 and 24 /i~^ Mpc. The richness for all r^i clusters is taken to be d = 90 h^^ Mpc. The heavy 
and thin lines are for OCDM and SCDM, respectively. 

Fig. 12. — The cross-correlation function between cluster and its particles for r^ = 6 h~^ Mpc 
clusters with richness d = 40 - 120 h^^ Mpc in the OCDM model. 

Fig. 13. — Temperature function n{T, r^i) averaged from 8 realizations for clusters of r^/ = 0.75, 
1.5, 3, 6 and 12 /j-^ Mpc in the SCDM, LCDM and OCDM models. The TFs of r^ = 12 /i^^Mpc 
(dashed lines) are clearly not a member of the scaling family. The observations of Henry & Arnaud 
(1991) and Edge et al. (1990) are illustrated as triangles and circles, respectively. 

Fig. 14. — Temperature functions averaged from 3 realizations for clusters of r^i = 0.75, 1.5, 3, 6 
and 12 h^^ Mpc in the SCDM model with COBE normalization {a% = 1.0). The observational 
points are the same as in Fig. 13. 
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